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Abstract. The gravitationally lensed quasar B1600+434 
(z = 1.61, wiy = 21.6) has been observed at the 2.56m 
Nordic Optical Telescope (NOT). In this Letter we report 
the discovery of an edge-on late-type galaxy located be- 
tween the two lensed components (separation l'/4) , close to 
the fainter image. The galaxy photometry indicates that 
its redshift is approximately 0.4. We detect a large colour 
difference between the two images due to significant ob- 
scuration of the faint image. The estimated amount of 
absorption as a function of colour indicates that the ex- 
tinction may be due to dust in the lensing galaxy. We also 
present evidence of flux variability in B1600+434 with a 
detected change of ^0.25 mag in one year. The theoreti- 
cally expected time delay is of the order of one month and 
so the system may be an interesting object for determining 
the Hubble constant. 
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1. Introduction 

As part of the CLASS radio survey Jackson et al. (1995) 
reported the discovery of a new gravitational lens (GL) 
candidate, B1600+434, with an image separation of 9 = 
l'/4 and a flux ratio of I A /I B = 1.3 ± 0.04 at 8.4 GHz. 
Traditionally, gravitationally lensed images are expected 
to have identical colours due to the achromatic nature of 
gravitational lensing. In this Letter we report our find- 
ings from optical multicolour observations of this system 
obtained at three epochs which reveal a large colour dif- 
ference between the two components. 

* Based on observations obtained with the Nordic Optical 
Telescope (NOT), La Palma 

** Present address: NORDITA, Blegdamsvej 17, DK-2100 
K0benhavn 0, Denmark (jens@nordita.dk) 



2. Observations and reductions 

The observing log is given in Table 1. For the observa- 
tions we used the CCD cameras BroCam 1 (epochs 1 and 
2) and BroCam 2 (epoch 3). BroCam 1 is equipped with a 
thinned backside illuminated TEK1024 CCD with a gain 
of q = 1.7e~ ADU -1 , a readout noise of RON = 6.5e~ and 
a pixel scale of 0'.'176 pixel -1 . BroCam 2 applied a thinned 
Loral 2048 CCD with q = 1.27c- ADU -1 , RON = 6.22e - 
and a pixel scale of O'.'ll pixel -1 . A finding chart of the 
object and stars in the field is shown in Figure 1, includ- 
ing the three stars used as internal calibration and point- 
spread function (PSF) stars. In what follows we label the 
Northern image of B1600+434 as image A and the fainter 
Southern image as image B (cf. Fig. 2). 



Table 1. Dates and combined exposure times 



Epoch 


Date 


B 


V 


R 


I 






sec 


sec 


sec 


sec 


1 


July 23-24 1995 


1200 


800 


600 


1200 


2 


September 10-12 1995 


1200 


2300 


600 


1200 


3 


June 7-10 1996 


1800 


2800 


900 


1200 



The CCD frames were overscan subtracted, zero- image 
subtracted and flat fielded. Each science frame was man- 
ually cleaned for cosmic rays in the vicinity of stars and 
objects used in the analysis. 

Standard stars in M5, M92 and a few Landolt fields 
were observed on June 7 1996 (epoch 3) under photomet- 
ric conditions and the three PSF stars could thereby be 
calibrated to the Johnson-Kron-Cousins system. This was 
achieved by using calibrated photometry from Sandquist 
et al. (1996) (M5), Davis (1995) (M92), Grundahl (1996) 
(M5, M92) and Landolt (1992). The transformation equa- 
tions were fitted to about 80 stars to derive colour and 
extinction terms. The calibrated standard star magni- 
tudes could be reproduced to an rms of 0.02 — 0.03 mag. 
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Fig. 1. Combined (73 x 83 arcsec) image of the B1600+434 
field. North is up and east is to the left 

The extinction coefficients determined at the Carlsberg- 
Automated-Meridian-Circle (CAMC) for the La Palma 
site are in good agreement with the values deduced from 
our observations. 

3. PSF photometry 

Using DAOPHOT II (Stetson 1987) we determined an an- 
alytical PSF with a residual lookup-table from the three 
bright stars in the field. The PSF photometry was carried 
out using ALLSTAR. A multiple fit of the PSF to the QSO 
pair was constrained by fixing the position of B relative 
to A using the radio positions. A background value simi- 
lar to the immediate surroundings was used. This yielded 
magnitudes of the QSO images relative to the PSF stars. 
The photometric zero-point was determined from aperture 
photometry of the PSF stars. The calibrated magnitudes 
for A, B and the PSF stars are given in Table 2 for each 
colour band in epoch 1, along with the average measured 
FWHM values of the three stars. 

4. Lensing galaxy 

After subtracting the scaled PSFs from image A and B 
we are able to study the residual lens galaxy located very 
close to image B (see Fig. 2). We estimate the lens to have 
a major axis a = 4'.'1 ± 0'.'2, an axis ratio a/b = 2.4 ± 0.2 
and a position angle PA = 46° ± 3°. These values were 



Table 2. Photometry (mag) of image A, B and stars (epoch 1) 





FWHM 


A 


B 


B-A 


1 


2 


3 


B 


1'.'21 


22.12 


23.82 


1.70 


18.97 


18.74 


18.73 


V 


l'.'Ol 


21.92 


23.32 


1.40 


17.69 


17.75 


17.87 


R 


0'.'60 


21.46 


22.34 


0.88 


16.81 


17.12 


17.29 


I 


0'.'77 


20.92 


21.43 


0.51 


16.09 


16.63 


16.84 



measured from an isophote 1.2cr above the sky using SEx- 
tractor (Bcrtin 1996). Astrometry of the lens, the South- 
East (SE) galaxy and QSO images relative to image A 
(RA 20 oo : 16 h 01 m 40!5, Dec 20 oo : +43°16'47'.'0) is given in 
Table 3. We apply an elongated aperture with a given ma- 
jor and minor axis to measure the intensity of the lensing 
galaxy. The sky is estimated from an outer elongated an- 
nulus extending from 1.5 to 3.0 times the major and minor 
axis. To prevent neighbouring galaxies from contaminat- 
ing the sky estimate, overlapping apertures of such objects 
were masked out. We also measured the two galaxies to the 
South-East (SE) and West (W) of B1600+434. To correct 
for neighbouring stars or galaxies surrounding PSF stars 
2 and 3 we measured these stars with a small aperture 
radius and added aperture corrections based on growth 
curve plots (Grundahl & Hjorth 1995) of star 1 (see Fig. 1) 
which had a nicely converging growth curve in all images. 
The photometric results are given in Table 4. 



Table 3. Astrometry of lens, SE galaxy and QSO images 





Aq 


AS 


A 


0.0 


0.0 


B 


Of! 72 


1'.'20 


Lens 


0'.'74±0'.'10 


0'.'85±0'.'10 


SE 


4" 33 ± Of! 05 


3727 ± 0705 



Table 4. Photometry (mag) of lensing galaxy 

B V R I 

Lens 23.6 22.0 21.1 20.3 

SE 22.0 20.6 19.9 19.2 

W 23.9 22.3 21.7 21.0 

Using the B — V, V — R and R — I colour indices we 
estimated the photometric redshifts of the three galaxies. 
These colours were fitted (by least squares) to the tem- 
plate colours of various redshifts and galaxy types tabu- 
lated in Fukugita et al. (1995). On the basis of the residual 
image and the measured galaxy properties we assume that 
the lensing galaxy is a spiral galaxy. The photometry indi- 
cates approximate redshifts of zi ~ 0.4 for the lens, 0.1 for 
the SE galaxy (E/S0) and 0.2 for the W galaxy (spiral). 
The estimated uncertainty in these redshifts is 0.1. 
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Fig. 2. Combined V, R and I image of B1600+434 showing the 
A and B components (left). The same field after image A and 
B have been subtracted reveals the morphology of the lensing 
galaxy more clearly (right). The early- type SE galaxy is also 
seen in the images. The frames measure 14 arcsec on each side. 
North is up and east is to the left 



5. Galaxy extinction/absorption 

A particularly interesting property of B1600+434 is the 
strong apparent colour difference between the two images 
(see Table 2). Assuming that the true intensity ratio of A 
and B is that measured in the radio (Jackson et al. 1995), 
we have 



I_A 

Ib 



I_A 

Ib 



= 1.3. 



(1) 



radio 



If in addition we assume that A is not affected by redden- 
ing, i.e., I Ax = I Atrue and I Bx = k^Is^, the degree of 
absorption in the galaxy is given by 



kx 



(2) 



at restframe wavelength A. In Table 5 we give the observed 
colour band, the central wavelength (CWL) bandpass, the 
restframe CWL assuming z = 0.4, the A:B intensity ra- 
tio and the estimated absorption values. The reddening of 



values do, however, roughly agree with the expected trend 
for dust obscuration of image B. Note that Nadau (1991) 
estimated the first differential reddening of a gravitational 
lens system (2237+030) and found good agreement with 
extinction due to dust in the Galaxy. 

6. Variability 

To measure the possible effects of variability of the QSO 
we have performed aperture photometry on the combined 
light of the lensed images (A+B), the lensing galaxy and 
the SE galaxy. This method was chosen because the see- 
ing was worse than one arcsec in several images making 
the PSF subtraction uncertain for variability studies. As 
a check on our results we also used an aperture covering 
only the A+B images and the lens and found similar re- 
sults, but with larger uncertainties (especially in the B 
band). To allow an optimal analysis of the three epochs, 
the images were magnified, shifted and rotated to the same 
coordinate system. A single circular aperture was centred 
between the SE and lens galaxies. Due to low signal-to- 
noise, errors in the sky estimate and contamination by 
neighbouring objects, the growth curve of the aperture 
did not converge at large radii. To overcome this problem 
we reduced the aperture radius to a minimum of 5'.'0, while 
still including most of the light. The flux outside this ra- 
dius was included by extrapolating the growth curves and 
determining an aperture correction for each image. The 
sky was estimated from a sky annulus with a radius of 
8'.'8 extending to 9'.'9. The measured total magnitude of 
A+B+the lens+the SE galaxy is listed relative to epoch 1 
in Table 6. The uncertainty in the B band magnitudes is 



Table 6. Variability of A+B+lens+SE in the B, V, R and I 
bands relative to epoch 1 





B 


V 


R 


I 


A 2 _i 


0.00 + 0.10 


0.11 + 0.02 


0.10 + 0.02 


0.18 + 0.02 


A 3 _i 


0.30 + 0.10 


0.25 ± 0.03 


0.15 + 0.02 


0.21 + 0.02 



Table 5. Lensing galaxy rest-frame absorption values 



Filter 


CWL 


CWL Z=0A 


Ia/Ib 


kx 


B 


4448 


3177 


4.79 


3.68 


V 


5505 


3932 


3.63 


2.79 


R 


6588 


4706 


2.25 


1.73 


I 


8060 


5757 


1.60 


1.23 



image B is most likely due to dust obscuration in the lens- 
ing galaxy. The extinction due to dust increases roughly as 
k\ oc 1/A (Mathis 1990). The R and I band k\ estimates 
follow this trend, while the B and V gradients are some- 
what higher than expected. This discrepancy is probably 
partly due to errors in our PSF photometry of the B and 
V colours because of the worse seeing in these images. The 



approximately 4-5 times greater than those in the other 
filters. Apart from this filter the remaining bands show a 
consistent variation. Due to the uncertainties in the mea- 
surements it is difficult to identify whether the detected 
variation is due to an intrinsic flux variation of the QSO or 
a micro-lensing event. Jackson et al. (1995) indicated flux 
variability in their radio data and this is clearly supported 
by our optical data. 

7. Lens mass and time-delay 

Borgeest & Refsdal (1984) introduced a radially symmet- 
ric galaxy model to give rough estimates for the mass and 
time-delay in the QSO 0957+561 system. We apply their 
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model to estimate the lensing mass and time-delay be- 
tween image A and B in B1600+434. The lensing galaxy 
is very close to the position of the B component and we 
found that the separation between image A or B and the 
lens along the vector connecting A and B is 9 a = 1" 1 1 and 
9b = 0'/29. The cosmological correction factor, T, is found 
in Kayser & Refsdal (1983) for the lens and source red- 
shifts and a flat Fricdmann universe model (qo = 0.5 and 
H a = 100ft km s -1 Mpc -1 is used throughout this Letter). 
We assume there is no additional shear from a cluster or 
other galaxies (7 = 0) and that the lens has an isother- 
mal mass distribution (k = 1). Inserting the observational 
values into equation (15b) of Borgeest & Refsdal (1984) 
gives a time-delay of 24 ft^ 1 days. The enclosed mass in- 
side the angular radius 9/2 (the Einstein radius) is found 
using their equation (11) and gives 6.8 ft -1 10 10 M©. 

As the lensing galaxy of this system is relatively bright 
it is straightforward to estimate the mass-to-light ratio. 
We measure the galaxy luminosity within a radius of 
9/2 (= 0'/7 = 1.9 kpc) from the galaxy center and com- 
pute the rest-frame absolute luminosity. The rest-frame 
luminosity of the galaxy is estimated using the method 
described in van Dokkum & Franx (1996); note that there 
is a sign error in their definition of the conversion con- 
stants Cj. From the observed V, R and I colours (given 
in Table 4) we find the rest-frame B and V magnitudes 
at z = 0.4 to be tub — 23.6 and my = 22.8. Assum- 
ing ft = 0.7 the corresponding absolute magnitudes are 
M B = -17.7 and M v = -18.6. This gives mass-to-light 
ratios of 51 and 43 in solar units for the rest-frame B and V 
bands, respectively. These relatively high M/L estimates 
are, however, affected by various effects. Edge-on spiral 
galaxies commonly have dust lanes which significantly ob- 
scure the luminous matter in the galaxy (Barnaby 1995). 
Furthermore, the seeing smears out the light inside the re- 
stricted aperture used and the PSF photometry unavoid- 
ably introduces errors with a tendency for oversubtraction. 
Our mass-to-light ratios should consequently be taken as 
upper limits. Infrared observations of B 1600+434 with 
HST would offer an optimal measurement of the true lens 
galaxy M/L. 

8. Conclusions 

We find that a V — 22.0 nearly edge-on spiral galaxy 
acts as the lens in B 1600+434 and causes a significant 
reddening of image B, due to dust obscuration. The red- 
dening of image B is unusually strong and may in general 
have implications on lensing statistics, because most opti- 
cally discovered GL candidates are chosen by their colour 
similarity (e.g. Jaunsen et al. 1995). Small deviations in 
the relative magnitude for a given filter is usually allowed, 
but a difference of the order of one magnitude or more will 
probably result in that object being discarded as a GL can- 
didate. The failure to detect all GL candidates in a sample 
will introduce systematic effects which affect the results of 



the statistical analysis. Kochanek (1991) discussed various 
bias effects in optical gravitational lens surveys and con- 
cluded that 10-30 % of all GL systems are lost already 
in the process of identifying the QSOs. However, this is a 
rough estimate which depends on the selection criteria of 
the QSO sample. Indeed, in view of the results presented 
here, the number of unidentified GL systems in optical 
QSO/GL surveys is difficult to estimate. 

If most of the GL systems discovered so far are pro- 
duced by early-type ellipticals this could be explained by 
the fact that optically thick spiral galaxies tend to redden 
one of the (two) lensed images, thereby causing many spi- 
ral lensed systems from being detected in optical surveys. 
Radio selected samples are much less affected by such se- 
lection effects and may consequently reveal more spiral 
lenses. In fact B1600+434 is the second radio-selected spi- 
ral lens galaxy discovered at NOT by our group, the first 
being B0218+357 (Grundahl & Hjorth 1995) and both 
systems show signs of significant extinction. 

Our variability analysis indicates a variation of about 
0.25 mag year -1 . The lens system B1600+434 may there- 
fore be an interesting candidate for determining the time- 
delay and Ho. We estimate the theoretically expected 
time-delay to be of the order of one month. It should be 
noted, though, that the shear effects of the SE galaxy may 
affect the lensing and should probably be accounted for in 
a realistic model of the lens potential. 
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